Confinement-deconfinement transition in a generalized Kitaev model. 
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We present a spin model, namely, the Kitaev model augmented by a loop term and perturbed by 
an Ising Hamiltonian and show that it exhibits both confinement-deconfinement transitions from 
spin liquid to antiferromagnetic/spin-chain/ferromagnetic phases and topological quantum phase 
transitions between gapped and gapless spin liquid phases. We develop a Fermionic mean-field 
theory to chart out the phase diagram of the model and estimate the stability of its spin liquid phases 
which might be relevant for attempts to realize the model in optical lattices. We also conjecture 
that some of the confinement-deconfinement transitions in the model, predicted to be first order 
within the mean-field theory, may become second order via a defect condensation mechanism. 
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PACS numbers: 74.45+c, 74.78.Na 

Quantum phase transition from ordered to paramag- 
netic phases in two-dimensional (2D) spin models has 
been a subject of recent interest [lj] . Such paramagnetic 
phases and associated quantum phase transitions have 
been conjectured to be of relevance to the properties 
of several strongly correlated systems including cuprates 
and quasi-2D organic materials 0, S, 0] • A class of these 
paramagnetic phases which do not break any constituent 
symmetries of their underlying lattice are called spin- 
liquids and are generally believed to be natural candi- 
dates for paramagnetic phases obtained by disordering 
non-collinear spin-ordered magnets [5j . However, the pre- 
cise criteria for realization of spin liquids and the nature 
of quantum phase transitions to them from spin-ordered 
phases are far from being settled issues (6j. The physics 
of these spin-liquid phases can be described by represent- 
ing the spins in terms of spinons, which are Fermionic 
CP(N) fields, coupled to bosonic gauge fields @, @. In 
most commonly studied examples, the symmetry group 
associated with these gauge fields are either U(l) or Zi\ 
the corresponding spin- liquids being dubbed as f/(l) or 
Z 2 spin liquids. In the ordered phase of spins, the spinons 
are confined while the spin-liquid paramagnets constitute 
phases with gapped or gapless deconfined spinon excita- 
tions. Thus the transition between these phases serve 
as examples of confinement-deconfinement transitions for 
spinons. Examples of such transitions has been stud- 
ied in several high energy and condensed matter models 

In this letter, we present a spin model, namely the 
Kitaev model augmented by a loop term and perturbed 
by an Ising model, on a 2D hexagonal lattice with the 
Hamiltonian 

H = H K + H L + Hr, H K = -J2J2 J «°i °& 

j£ A a link 

H * = AJ E E h l = -kJ^2w p (l) 

j all links p 



where j a is the nearest neighbor of j connected by the 
a = x,y, z link of the hexagonal lattice as shown in Fig. 
[TJ a a are the usual Pauli matrices, W p — <rf afcrf <rf <rf <t| 
is the loop operator, where, as shown in Fig. [1] 1..6 de- 
notes sites of a hexagonal plaquette, A and B denotes two 
sublattices of the honeycomb lattice, and j' denotes all 
nearest neighbors of j. We demonstrate that this model 
exhibits confinement-deconfinement transitions from de- 
confined spin liquid to confined Ising ordered antiferro- 
magnetic (AFM), ferromagnetic (FM) or spin-chain (SC) 
phases and can therefore serve as a test bed to study 
such transitions. In addition, we show that the model 
also supports two distinct spin liquid phases with gapped 
and gapless deconfined spinon excitations and exhibits a 
topological quantum phase transition between them. To 
the best of our knowledge, the model presented here is 
the only spin model which supports spin liquid phases 
with both gapped and gapless deconfined spinon excita- 
tions and displays both confinement-deconfinement and 
topological quantum phase transitions. We chart out the 
phase diagram of this model using a Fermionic mean- 
field theory and estimate the stability of the deconfined 
phases as a function of the strengths of the loop and 
Ising terms. There have been proposals for experimen- 
tally realizing the Kitaev model in systems of ultracold 
atoms and molecules trapped in optical lattices [HI . Any 
such physical realization of this model will always have 
contaminating interactions; our work thus provides an 
estimate of the stability of its gapless spin liquid phase 
against such interactions. Finally, we conjecture that 
the confinement-deconfinement transition between the 
gapped deconfined Kitaev phase and the confined Ising 
AFM or FM phases, predicted to be first order within 
the mean-field theory, may turn out to be second order 
via a defect condensation mechanism at large k. 

The Kitaev model with Hamiltonian Hk is a rare ex- 
ample of a 2D spi n model which can be exactly solved 
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FIG. 1: (Color online) Schematic representation of the Kitaev 
model on a honeycomb lattice showings the different links x, 
y and z and the two sublattices A and B. The sites labeled 
1..6 and their spin configuration is a schematic representation 
for a classical configuration with {W p ) — 1 



the parameter regime L/i — J2I < ^3 < Ji + J2, sup- 
ports a gapless phase [10( which can be described as a Z2 
spin liquid constituting a Fermi sea of deconfined gapless 
spinons and static Z 2 gauge fields. These spinons be- 
come gapped beyond this regime. Both the gapped and 
the gapless spinon phases have unit expectation value for 
a non-local loop operator W p : (W p ) = 1. They are dis- 
tinguished by the value of the Chern number which is ±1 
for the gapless phase and for the gapped phases. 

In the present work, we study the effect of augment- 
ing the Kitaev model with a loop term and perturbing 
it with an Ising Hamiltonian. In what follows, we shall 
study the Kitaev model in the isotropic limit and scale 
all energies by J: J a = J = 1 and obtain the phase dia- 
gram of H as functions of dimensionless couplings k and 
A. Wc note at the outset, that the ground state of the 
system at A — is the Kitaev ground state with gapless 
deconfined Fermionic spinons coupled to a frozen gauge 
field configuration. This ground state is exactly captured 
by the Fermionic theory which makes it an ideal starting 
point for our analysis. The ground state of the system 
at A — > ±00 is the Ising AFM/FM which are magnet- 
ically ordered phases with confined spinons. Thus the 
model, by construction, clearly supports confmement- 
deconfincment transitions. 

We begin by mapping the spin model H to a Fermionic 
model Hp, by using the Jordan- Wigner transformations 
[13]. We choose a path {i(n)} which runs along the 
x and v bonds and define a z , > = icA/„\b z ., .,, a x , n \ = 

C j(2n) rim<2n CT fe(m)> Cr j(2n+1) = ^j(2n+l) IIm<2ri+l G k{m) ' 

and a j = iv'jPj, where Cj, 6| are Major ana fermions 
operators at site j [(j(2n),j(2n — 1)) are the x bonds 
and (j(2n + 1), j(2n)) are the y bonds]. The resultant 



Hamiltonian becomes 

Hp = -J2 [ E ic J c j' a + E ''vv'y'v 

jGA OL—x.y links z link 

j,k£plaquette zlink 

+ A E E ''yy'v </ ( 2 ) 

j a— all links 

where the subscript j, k G plaquette indicates that the 
sum is over sites which belong to the A sublattice of a 
given plaquette as schematically shown in Fig. [TJ Note 
that for A = 0, the operators ibjb*, , commute with the 
Hamiltonian and are therefore a constants of motion. In 
this limit, H is exactly solvable. When A is turned on, 
these operators acquire dynamics and their fluctuations 
arc ultimately expected to confine the spinons through a 
confincmcnt-dcconfincmcnt transition. 

To make further progress, we treat Hp within an RVB 
type mean-field theory [3] and introduce the mean-fields 
on the sites (corresponding to spin ordering) and on 
links (corresponding to the emergent gauge fields) of the 
hexagonal lattice: (ibjCj) = (crj) = A 1 ( 2 ), (ibijCf ) = Pa, 
(ibjb z , ) = 7 Q , and (ic,-c-/ ) = 7o Q . Note that keeping in 
mind the bipartite nature of the hexagonal lattice and to 
allow for possible AFM phases, we have introduced two 
mean-fields Ai and A2 corresponding to the two sublat- 
tices shown in Fig. [TJ In terms of these mean-fields, one 
can decompose the quartic terms in Hp as 

ibplibib* K = i^bl+bibl,)-^ 

ibjCjib^Cj^ = JaiCjC^ + jQaibjb*, - Ja'yoa 

-Aiife^c^ - A 2 ib z j c j + AiA 2 
-iMb^+b^-Pl (3) 

where we have only considered mean-fields which are ei- 
ther on-site or reside on links between two neighboring 
sites. The resultant quadratic mean-field Hamiltonian, 
in momentum space, can be written as 

k 

+J (p - 2k 1z /J + e lkl + e z( - kl+k2 ^ bfb% 

+(icf6f -i&fcf)(&(l + A)+&e*i 

+/3 y e lk -) - Cl bfc£ - c 2 6f f cf + h.c 

+k 7z 2 - (1 + X) lzl0z + (1 + 3A)AxA 2 (4) 

where J = (1 + Xjx), aJ Q = (1 + A)7 Z , j' = \~fQ X ), 
/3J = (1 + A) 7oz, c i(2) = (1 + 3A)A 1 ( 2 ), and the momen- 
tum k = k\e\ + ki&i with the unit vectors ei = x + y/V3 
and &2 — 2yj\f%. 

We now minimize H m { numerically and obtain the 
mean-field phase diagram of the model as a function of 
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FIG. 2: (Color online) The mean-field phase diagram for 
the model. The blue dashed lines represent confinement- 
deconfinement transitions which may become second order 
via defect condensation mechanism. The triple points occur 
at (A*,«*) = (1.01,1.07) and (A^kS) = (3,5.64) for A > 
and (A3, K3) = (-1.01, 1.88) for A < 0. 



A and k. Note that the mean field solution is exact at 
A = 0. This phase diagram is shown in Fig. [2j In ac- 
cordance with our earlier expectation, we find that at 
large positive (negative) A, the ground state of the model 
is an Ising AFM (FM) which corresponds to confined 
phase of spinons while at small A, the model exhibits a 
deconfined gapless phase DC1. The transition between 
these two phases at low k and for positive A occurs via a 
SC phase, which corresponds to AFM alignment of spins 
along chains in x direction of the hexagonal lattice with 
ferromagnetic arrangement of such chains in the y di- 
rection. For negative A, there is a direct transition to 
the FM phase. At high enough values of n, we find 
another gapped deconfined phase DC2 and second or- 
der topological quantum phase transitions between DC1 
and DC2 phases. The confinement-deconfinement tran- 
sitions at high k always occur from DC2 to AFM/FM 
phases. These transitions are predicted to be first or- 
der within mean-field theory. The phase diagram ex- 
hibits two triple points at (A*,«*) = (1.01,1.07) and 
(A|,k|) = (3,5.64) for A > 0. These represent meeting 
points of AFM, SC and DC1 and AFM, DC2 and DC1 
phases respectively. For A < 0, there is one triple point 
(A|,/t|) = (-1.01,1.88) where the FM, DC1, and DC2 
phases meet. We also note that our mean-field analysis 
also gives an estimate for the stability of the deconfined 
phase of the Kitaev model (-0.07 < A c < 0.08 for k = 0) 
under external perturbing Ising term which may be im- 
portant for physical realization of the Kitaev model and 
for quantum computing proposals based on it (l5. 17|. 

The plot of the loop order parameter (W p ), the spinon 
gap, and the AFM and the FM order parameters as ob- 
tained from the mean-field theory, is shown, for k = 7, 
as a function of A in Fig. [3] We note that all the order 
parameters show discontinuous changes at the transition 
points indicating first order transitions. The spinon gap, 
in contrast, increases linearly and continually with A in- 
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FIG. 3: (Color online) Plot of the loop order parameter (solid 
blue line), the spinon gap (red dashed line) and the FM and 
the AFM order parameters (black dash-dotted lines) as a 
function of A for n — 7. 



dicating a second order quantum phase transition be- 
tween DC1 and DC2 phases. The presence of this topo- 
logical quantum phase transition and the linear varia- 
tion of the spinon gap with A can be understood qual- 
itatively from H m {. For large k, it requires a large A 
to destabilize the Kitaev ground state in favor of Ising 
AFM/FM. In addition, numerically we find that in the 
Kitaev phase 72(7^) ~ 1(0). As a result, beyond a crit- 
ical value of A = A c , the effective couplings along the 
links, J1.2 ~ (1 + A7 X ), J3 ~ 7z(l + A), fail to satisfy 
I J\ — J2 1 < J3 < J\ + J2 thus leading to a gapp ed phase 
via a topological quantum phase transition [lOl [Tl| . The 
spinon gap in this gapped phase varies linearly with J3 
[ill [Til l and hence shows a linear variation on A. At 
small k, the confinement-deconfinement transitions to 
the SC/FM phases occur before A c is reached and hence 
the topological phase transitions do not occur. 

Next, we argue that it is possible that some of the 
confinement-deconfinement transitions, predicted to be 
first order by mean-field theory, can become second order 
in the large k limit via a defect condensation mechanism 
[i~8] |. We note that similar second-order transitions have 
been seen in numerical studies of related models 1(| Tj^ . 
We begin with the case A > and start from the Ising 
AFM ground state. The possible classical defect spin 
configurations over this ground states whose condensa- 
tion can lead to either the Kitaev or the SC phases are 
shown in Fig. [4] The energy of these defects are given 
by E sc = 8(A - 1), E DC i = 14A - 6 - k, and E DC2 = 
10A-2-K. Note that E AFM = E DC i = E sc = E DC 2 = 
at A = A* = 1 and k = k* — 8 which represents a triple 
point in the n — A plane. For the SC phase to occur, we 
need E sc = and E sc < Edci, Edc2 which yields the 
conditions 



X s c =1 < c < 6A + 2, nf < 2A + 6 (5) 



Similar analysis lead to conditions for instability due to 
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FIG. 4: Schematic representation of classical spin configura- 
tion of the AFM phase and possible defects that can destabi- 
lize this phase. See text for details. 



condensation of the defects DC1 and DC2 
k° c1 3 



A; 



DO 
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,DC2 
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X? C1 < 1, 



A 



DC2 



> 1, 



DO 



„DC2 



> 6A 



> 2A 



2,(6) 
6.(7) 



From Eqs.[5l[6j and[71 we find that the SC phase desta- 
bilize the AFM phase at A = 1 for all k < n* . This 
represents an order-order transition which is expected to 
be first order. As we decrease A further below 1, the de- 
confined gapless Kitaev phase takes over at some point. 
The nature of this transition is not predictable from this 
condensation mechanism, since the transition do not take 
place from the AFM phase. For k> k* , the AFM phase 
is destabilized in favor of the deconfined phase DC2 when 
A c = k c /10+1/5. The DC2 phase represents a disordered 
phase in terms of the spin variables and thus this transi- 
tion can be second order. For k > «*, the energy of DC1 
becomes lower than that of DC2 for A < 1. This indi- 
cates a second order transition between the two decon- 
fined phases with no broken symmetries and is therefore 
a the topological quantum phase transition. Note that 
this transition line is independent of k and therefore is 
expected to appear as a vertical line in the n — A plane. A 
similar discussion for A < shows that the defect conden- 
sation energies of the phases DC1 and DC2 (the SC phase 
is not energetically favorable for A < 0) over the FM state 



are E' T 



DO 



14|A| + 6-k and E' 



" DC2 = 10|A| + 2-k. We 
find that such a mechanism works only for large k where 
E' DC1 and E' Dc2 can be negative and in this regime the 
FM phase is always destabilized by the DC2 phase lead- 
ing to a confinement-deconfinement at |A| = k/10 — 1/5. 
For lower k, we expect to find a transition between the 
DC1 and the FM phase which is not mediated by such 
defect condensation. The phase diagram obtained by this 
simple qualitative discussion of the defect condensation 
mechanism described above has most of the major fea- 
tures shown in Fig. [5] and a direct AFM to DC1 transi- 
tion line at intermediate k for A > 0) and this gives us 



some confidence about it's qualitative correctness [20j. A 
quantitatively accurate study of the obtained phase dia- 
gram requires a more sophisticated treatment of quantum 
fluctuations and is left for future study. 

To conclude, we have presented a spin model, namely 
the Kitaev model, augmented by a loop term and 
perturbed by an Ising Hamiltonian, and have shown 
that the model exhibits a rich phase diagram with 
both confinement-deconfinement and topological quan- 
tum phase transitions. We have estimated that the topo- 
logical phase of the Kitaev model is unstable to about 
10% contamination by Ising interactions. 
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